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Abstract We calculate the Casimir force at a finite cut-off A by summing the forces induced by the all fluctuation 
modes. We show that the Casimir force is independent of the cut-off function in the limit LA — oo. There is a correction 
in the order of (LA)~?, when LA is finite and large. This correction becomes remarkable when L is comparable with 
the microscopic length scale A~'. It has been demonstrated that the Casimir force at a finite cut-off should be defined 
by summing forces of all fluctuation modes, instead of the derivative of Casimir energy with respect to L where an 
additional derivative of the cut-off function has been introduced. 
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1 Introduction 

In 1948, Casimir predicted an attractive force between 
parallel conducting plates due to vacuum fluctuations of 
electromagnetic (EM) field."] The study of Casimir effect 
has been a subject of extensive research, especially after 
significant progress in experiments?—*) (for a review on 
Casimir effect, see Ref. [6]). 

Because the vacuum energy is notoriously ultraviolet- 
divergent, one must use some regularization methods to 
extract the finite Casimir energy and define the Casimir 
force as the derivative of Casimir energy with respect to 
system size L. The most useful regularization methods 
are the zeta functional regularization method!) and the 
cut-off technique. 

The zeta functional regularization has a pure mathe- 
matical character. The cut-off technique is of a physical 
character. With finite cut-off all the physical quantities 
The dropped parts of the vac- 
uum energy, which are divergent when A — ov, are the 


are free of divergence. 


bulk, hyper-surface,---, surface and the edge energies’! in 
a rectangular box. It has been shown that the Casimir 
energy is independent of the cut-off function when the 
cut-off A — oo.8-1° Secondly, for quasi-one-dimensional 
systems such as film and piston system, it has been found 
out that the divergent parts inside and outside cancel each 
other and do not contribute to the Casimir force.{!!~15] 
Physically there should be no real divergence. Some 
kinds of cut-off functions must exit at ultra-high energy 
scale. In principle, quantum field theory (QFT) should be 
regarded as a low-energy effective theory of a more fun- 
damental theory. When we use the QFT to the ultra- 
high energy scale, as we do in calculating the vacuum 
energy, we should introduce a cut-off function C(|k|/A) 
to describe the modification of QFT at high energy scale 


in principle. As an example, one can take a sharp cut- 
off related to Plank energy scale because of the smallest 
scale in the space-time. Another example is the Gaus- 
sian damping function as a result of the space-time non- 
commutativity.!!6] Of course, the exact form of the cut-off 
function depends on the unknown fundamental theory at 
high energy scale. 

Recently the Casimir energy at a finite cut-off was 
studied in Refs. [16] and [17]. The effect of the finite cut- 
off cannot be neglected when system is not large enough. 
For example, the correction is about 5% when L = 10471 
(Ref. [16]). For small enough system, the Casimir energy 
deviates from the standard behavior L74 and have a min- 
imum around L = A7!. By the derivative of the Casimir 
energy they define the Casimir force which flips its sign 
at the point of minimum.!!7] In the Casimir force they de- 
fined there is an artificial force related to the derivative 
of cut-off function C(|k|/A) with respect to L. At a finite 
cut-off we can obtain the bulk force of the vacuum. The 
derivative of the cut-off function gives also an additional 
contribution to the bulk force so that it becomes attractive 
and makes the vacuum unstable. This is also artificial. 

In this paper we calculate the vacuum force by sum- 
ming directly the forces induced by all vacuum fluctuation 
modes with a finite cut-off. This method can give a posi- 
tive bulk force, which makes the vacuum stable. From the 
vacuum forces inside and outside the wall we can get the 
Casimir force at a finite cut-off, which is different from the 
derivative of the Casimir energy. 

Our paper is organized as follows. In Sec. 2 we dis- 
cuss the two definitions of the Casimir force we mentioned 
above. In Sec. 3 we investigate the asymptotic behavior of 
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the Casimir force for large system size L >> AT! with gen- 
eral cut-off function. In Sec. 4, with exponential (4.1) and 
power (4.2) damping cut-off function the Casimir force is 
calculated for general L. Finally we make some conclu- 
sions and discussions in Sec. 5. 


2 Casimir Force at Finite Cut-off 

Without loss of generality, we study a massless 
scalar field y(t,2) confined in a one-dimensional (1D) 
interval 0 < a2 < L. Here we only consider typ- 
ical Dirichlet-Dirichlet (DD) boundary conditions with 
p(t,0) = y(t,L) = 0 and Dirichlet-Neumann (DN) 
boundary conditions with y(t,0) = O,p(t,x)|z=1 = 0. 
It is easy to obtain the allowable fluctuation modes 


kin(L) = = (n+ 8), (n =0,1,2,...), (1) 


where Bopp = 1,8pn = 1/2. The other idealized bound- 
ary conditions (such as Neumann—Neumann (NN), period- 
ical (PR), and anti-periodical (AP) boundary conditions) 
are quite similar to DD or DN boundary conditions.|®] 

The energy of a vacuum fluctuation mode k is given 
by (A=c=1) 

solk) = tolk) = E. (2) 

The force induced by this fluctuation mode can be ob- 
tained by the derivative, 


folLyk) = -Zeol = É, (3) 
which can be considered also as the radiation pressure of 
a virtual particle with momentum k. 

In literatures, !!-68-9:16.17] the Casimir force was defined 
mostly by the derivative of Casimir energy with respect to 
L. At a finite cut-off this definition gives an artificial force 
related to the the derivative of the cut-off function. So the 
Casimir force should be calculated by summing all forces 
induced by vacuum fluctuation modes. In this section the 
difference between these two definitions will be discussed. 


2.1 Casimir Force Defined by Summing Forces of 
Fluctuation Modes 


After summing the force fo(Z,k) of all fluctuation 
modes, we can get the vacuum force inside the interval. 
There is a notorious divergence if the summation is done 
until to infinite high energy scale. This divergence is arti- 
ficial because the quantum field theory used above cannot 
be applied to infinite high energy and should be modi- 
fied at very high energy scale. To take account of this 
modification, we introduce a dimensionless cut-off func- 
tion C(\k|/A) with 


C(0)=1, C(t 0) ~ o(x7’), (4) 


where C(0) = 1 means that there is no modification at 
low energy scale and C(x) decays faster than x7? so that 
the total force is convergent. 
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Then the vacuum force inside the interval is 
X O€o(kn) ~/kn\ 1% ka 
== kin A dD 
S o(a Enol), © 


called the modes-summing vacuum force 


Fo(L, A)= 


which is 
(MSVF). 

After taking the bulk limit L — oo, we obtain the 
modes-summing bulk vacuum force 


P(A) = lim Fo(L,A) = RA’, (6) 
where 
= °° dk 
R= Soe). (7) 
0 2T 


The fluctuation modes outside the interval gives the 
the bulk vacuum force. The Casimir force is obtained by 
the sum of the vacuum force inside and outside the interval 
kn 


w) AA. (8) 


F.(L,A) = m- F, = x ` knC( 
n=0 


We can rewrite the Casimir force in a scaling form 


F,(L, A) = L? F(A); (9) 
with the scaling function 
~ OT = n+ eas 
RM = 7 Lint Aol Te Fe. (10) 


2.2 Casimir Force Defined by Derivative of 
Vacuum Energy 


As a comparison, we discuss here the Casimir force de- 
fined by the derivative of the vacuum energy. At first the 
vacuum force inside the interval is 


aOR [kn] 
Compared with Eq. (5), this vacuum force has an addi- 
tional part, 


ee ee ey T 


where C() (2) = dC(a)/da is the derivative of cut-off 
function. The additional force is resulted from the deriva- 
tive of the cut-off function and should be artificial. 

We can argue this artificiality in analogy to the statisti- 
cal physics. If we take the cut-off function as a probability 
function, the present system is an analogy to a thermody- 
namic system. In thermodynamics the pressure is defined 
by the derivative of free energy rather than the internal 
energy. The derivative of free energy is equivalent to sum- 
ming forces induced by all particles. So the vacuum force 
at a finite cut-off should be calculated by summing forces 
of all fluctuation modes. 
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The bulk limit of the additional part of the vacuum 
force defined by the derivative is 


SF, (A) = lim 6Fo(L, A) = —2F A? . (13) 


The bulk vacuum force F, is positive and repulsive. But 
the bulk vacuum force defined by the derivative 


P(A) = P(A) + 8P (A) = -7 (A), 


which is negative and then attractive. The attractive bulk 
vacuum force makes the vacuum unstable and crash to- 
gether. This is an artificial effect. 

After subtracting the bulk vacuum fluctuation force we 
get the derivative Casimir force (DCF), 


(14) 


F.(L, A) = Fo(L, A) = Fy(L, A) 
i k k k = 
— n n (1) n 2 
TE En Pl) re (F) +AA 08) 
It can be rewritten also in a scaling form, 
F.(L, A) = L? EN), (16) 
where the scaling function 
F.(l) = F(t) + 6F.(0) (17) 
with an artificial part 
= L€ n+ = 
F) = = Ka yo eam GE 
dF.) = Un + AC (Te) t2. (18) 


In the following we will calculate the scaling function 
F.(1) of the mode-summing Casimir force. The artificial 
part ôF, (L) of the derivative Casimir force will be discussed 
also. 


3 Asymptotic Behavior of the Casimir Force 
for General Cut-Off 


Using the Able—-Plana sum formula, 


oA g g(0 
Yolo (n + 8)) =f dyg(y) + Gp. 
Qa Jo 2 
oa glay) g(—iay) 
+i | dy| e2m(yti8) 1 e2my—iB) — =k (19) 
we can rewrite the scaling function of Eq. (10) as 
7 C(it/1) C(-it/1) 
= o=- att] l= (t+iBr) _ 1 e2(t-iBr) — | -(20) 


From this expression we can obtain the asymptotic be- 
havior of this scaling function without knowing the specific 
form of the cut-off function C(x). 

3.1 Limit LA — oo 
Taking the limit 1 = LA — œ in Eq. (20), we have 


C(+it/l) = C(0) = 1 and obtain 

F.(00) = -7 Ba(8), (21) 
where B;(x) is the Bernoulli polynomial with 

EC ee ee ae (22) 
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Now the Casimir force can be written as 


F,(L, o0) = L~?F,(00) = —7 Bo) L~? 2 8 


This result agrees with the Casimir force calculated by the 
derivative of vacuum energy with cut-off A — oo. 

The sign of this Casimir force depends on the param- 
eter 8. For the Dirichlet—Dirichlet boundary condition 


FPD(L, o0) = eo <0. (24) 


The Casimir force in the DD boundary condition is at- 
tractive. For the Dirichlet-Neumann boundary condition 
FPN (L, œ) = ZL? > 0. 


1 1 
aG) =- 
2 12 48 


Therefore the Casimir force in the DN boundary condition 


(25) 


is repulsive. 


3.2 Large LA® 1 

For finite 1 = LA the scaling function F,(1) has correc- 
tion to its limit F.(co). Here we discuss the asymptotic 
behavior of this correction for l > 1. 


Because of the damping factor (e? + 1)~! in the in- 
tegrand of Eq. (20), for l > 1 we can expand C(+it/I) 
Correspondingly the scaling function be- 


around zero. 


comes 
F.(l) = Felo) + XO Foal, (26) 
i=1 
where 
Fai = — aC Oj Bagal), (27) 
il? ea) | E 
with the n-th derivative of the cut-off function 
d”C(x) 
c™ (0) = ; 28 
(0) T Maes (28) 


In Eq. (26) there are only even powers of I~! because 
3 = 1,1/2 for DD and DN boundary conditions. For gen- 
eral boundary conditions with other @ there should be also 
terms with odd powers. 

The coefficient of the leading correction is 


3 
Fog = — 7,0 (0) Bal), (29) 


which depends on the cut-off function (the second deriva- 
tive of it). 

When setting C®)(0) = 1, 1 = 5, the correction is 
about —1.97% for DD boundary conditions and about 
—3.45% for DN boundary conditions. C?)(0) could be 
larger than 1. For the power damping cut-off function 
used in subsection 4.2, C(?)(0) is of the order 10. Then 
the correction is remarkable for L ~ 5471. 

When the system size L becomes comparable with the 
microscopic length scale AT}, the higher-order terms of 
the correction must be considered. Then the correction 
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depends strongly on the form of the cut-off function. Al- 
though the condition L ~ A`! is hard to be achieved 
in the laboratory currently, it is possible for some theo- 
retical models, such as the Casimir’s semi-classic electron 
model!!®) and the bag model of quark.!!9] 


3.3 Large LA Behavior of Derivative Casimir 
Force 


We have got the general asymptotic behavior of the 
Casimir force calculated by summing forces of all fluctu- 
ation modes. Here we give the large LA behavior of the 
derivative Casimir force. 

The scaling function of the derivative Casimir force 
can be expanded similarly as 


Fe(l) = Felo) + X Feil, (30) 
j=l 
where 
Fosi = (1 T 2i) Fe oi . (31) 


When calculating the Casimir force, the artificial bulk 
vacuum forces inside and outside the interval cancel each 
other exactly. In the limit | = LA — oo, F.(l) gives the 
same result as we obtained in the last subsection. 

For large 1 = LA > 1 the correction terms of the 
derivative Casimir force are obviously larger than the cor- 
responding terms of the mode-summing Casimir force. 


4 Casimir Force for Specific Cut-off Functions 
To obtain the scaling function of the Casimir force in 
whole range, we need to know the specific form of the 
cut-off function, which depends on the physical properties 
at ultra-high energy scale. For example, the cut-off func- 
tion is a Gaussian damping factor when the space-time is 
noncommutativel!® and is quite similar in the presence of 
a minimal length.|° 
Here we consider an exponential damping cut-off func- 


tion 
C.(a4) =e”, (32) 
and a power damping cut-off function 
Cala) = (1+2), (s>2). (33) 


4.1 Exponential Damping Cut-off Function 


For the exponential damping cut-off function Ce(x), 
we obtain the scaling function of the Casimir force 
~ T | e7(1-6)/l ( 1 


#0 =F e eaa a |e 


For the DD boundary condition 6 = 1 and the scaling 
function becomes 


reo = Fam E) - (3) 


Because sinh(x) > x for any x > 0, so FP? (I) is always 
negative for any system size. When l > 1, the scaling 


(35) 
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function can be expanded as 
2 
FPP) = -2i Og |: 
Pey =-2ji- Ti+]. (36) 


The leading correction term here is consistent with the 
general expansion (27). 
In Fig. 1 we plot the correction RPP (1) defined by 


FPP) = FP? (oo) [1 + R?? (0). (37) 


We can see that the cut-off dependent correction is re- 
markable (25%) when L<3A7?. 


0.0 


—0.2 


—0.4 


R(LA) 


LA 


Fig. 1 Correction of the Casimir force in the Dirichlet- 
Dirichlet boundary condition and with an exponential 
damping cut-off function C.(|k|/A) = e7!*!/4, which is 
given in Eq. (37). 


We can get also the scaling function of the derivative 
Casimir force 


2 
FPP (i) = Z sinh? (Z) [1 F eoth( Z )| H L : 
This expression agrees with Eq. (6) of Ref. [17]. We find 
that FP? (1) does change its sign at l ~ 1, which has been 
claimed in Ref. [17]. But this change of sign is an artifi- 
cial effect of the derivative Casimir force. So the derivative 
vacuum force has artificiality not only in bulk, but also for 
finite LA. 
For the Dirichlet-Neumann boundary condition with 
(Bpn = 1/2), we get the scaling function of the Casimir 


(38) 


force, 
DN (1) — © gimh-2( © my Ë 
FPN (1) = Z sink (+) coth(=) = 689 
which can be expanded as 
~ Tr? 
FPN) =Z roa 4 
Psl- Br +oc4], 0 


for 1 >> 1. 
In Fig. 2 we plot the correction RPN (1) defined by 


FON (l) = FPN (oo)[1+ RPN (D). (41) 


Similar to the DD boundary condition, the finite cut-off 
correction is remarkable (25%) when L<4A7}. 
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With the exponential cut-off function we have shown 
that the finite cut-off must be considered when L is com- 
parable with the microscopic length scale A~!. In the next 
subsection we will consider the finite cut-off effects of the 
Casimir force with power damping cut-off function. 


0.0 
—0.2 


—0.4 


R(LA) 


—0.6 


—0.8 


—1.0 


Fig. 2 Correction of the Casimir force in the DN bound- 
ary condition and with an exponential damping cut-off 
function C.(|k|/A) = e7!*!/4, which is given in Eq. (41). 


4.2 Power Damping Cut-Off 


From the power damping cut-off function Cp(x), which 
has been given in Eq. (33), we obtain the scaling function, 


Fat) = I(E) ele- 4) -efe 


Z (s? — 38 4 uO (42) 

where E 
C(s, a) = X (nta) (43) 

n=0 


is the generalized Zeta function. 

The power damping cut-off function C,(a) is controlled 
explicitly by the parameter s. So it is convenient for us 
to study dependence of the Casimir force on the cut-off 
function. 

For the Dirichlet—Dirichlet boundary condition with 
Bpp = 1, the scaling function can be expanded as 


2 


T T j Z2 —4 
zi gt + WI? +O). 


The result of FP? (co) here agrees with the previous re- 
sults. Only the corrections of the Casimir force depends 
on the cut-off function, as we have discussed in Sec. 3 for 
general cut-off functions. 

We plot the correction of the Casimir force for the 
power damping cut-off function with s = 3,4,5 in Fig. 3. 
We can see that the Casimir force is not sensitive to the 
parameter s when l > 20, while dependent strongly on s 
when / < 10. The correction is larger for a faster damping 
cut-off function with a larger s. In comparison with the 


FPP (I) = (44) 


exponential cut-off function, the correction here is larger. 
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This can be explained by the larger second derivative of 
cut-off function Co (0) = s(s + 1), which is much larger 
than CP (0) = 1. 


0.0 


—0.2 


—0.4 


R(LA) 


—0.6 | 


—0.8 


1.0 i a 1 . fi 1 1 
0 5 10 15 20 25 
LA 


Fig. 3 Correction of the Casimir force in DD bound- 
ary condition and with a power damping cut-off func- 
tion Cp(|k|/A) = (1 + |k|/A)~°, which is obtained from 
Eq. (42). 


Similarly the scaling function for the Dirichlet- 
Neumann boundary condition can be expanded for l >> 1 
as 

FON (0) 


Dm | Tr? (45) 


Z 9 a +1)? + OT$). 
In Fig. 4 we plot the correction of the Casimir force in 
DN boundary condition for s = 3,4,5. Similar to the DD 
boundary condition, the Casimir force here is dependent 
on the cut-off function for small system, while is unsensi- 


tive to the form of the cut-off function for large system. 


0.0 


—0.2 
—0.4 | 


—0.6 


R(LA) 


—0.8 


1.0 bf 


0 5 10 15 20 25 
LA 


Fig. 4 Correction of the Casimir force under DN bound- 
ary condition and with the power damping cut-off func- 
tion Cp(|k|/A) = (1 + |k|/A)~°, which is obtained from 
Eq. (42). 


5 Conclusions and Discussions 

In the presence of a finite cut-off A, we have shown 
that the Casimir force should be defined by summing the 
forces induced by all vacuum fluctuation modes. The 
usual definition of the Casimir force by the derivative of 
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the vacuum energy has an attractive bulk force, which 
makes the vacuum unstable and is therefore unphysical. 
We have calculated the 1-dimensional Casimir force with 
the Dirichlet—Dirichlet and Dirichlet-Neumann boundary 
conditions. Its asymptotic behavior at large system size 
has been studied for general cut-off function. The com- 
plete behavior of the Casimir force for exponential and 
power damping cut-off functions have been investigated 
also. 

In limit LA — oo, it is shown that the Casimir force 
is independent of the cut-off function and in agreement 
with the result obtained by the derivative method. For 
large LA, there is the correction with a leading term pro- 
portional to (LA)~?. The derivative method has a strong 
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overestimation of the correction. 

The correction increases with the decrease of the sys- 
tem size L. When L is comparable with the microscopic 
scale A~!, the Casimir force deviates from the standard 
power law behavior and depends strongly on the cut-off 
function. Although the condition L ~ AT! cannot be 
achieved normally in laboratory, it could be possible for 
physical systems in microscopic scale such as Casimir’s 
semiclassical electron model.!!8] In the theoretical models 
of elementary particles the finite cut-off effect should be 
very important also. 

For the Casimir force of general d-dimensional systems, 
we expect that there are the similar cut-off dependent cor- 
rections for finite LA. 
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